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The reaction pp — > AA — > p-K + pn~ has been measured with high statistics at a beam momentum 
of pp — 1.637 GeV/c. The use of a transversely-polarized frozen-spin target combined with the 
self-analyzing property of A/ A decay allows access to unprecedented information on the spin struc- 
ture of the interaction. The most general spin-scattering matrix can be written in terms of eleven 
real parameters for each bin of scattering angle, each of these parameters is determined with rea- 
sonable precision. From these results all conceivable spin-correlations are determined with inherent 
self-consistency. Good agreement is found with the few previously existing measurements of spin 
observables in pp —* AA near this energy. Existing theoretical models do not give good predictions 
for those spin-observables that had not been previously measured. 

PACS numbers: 24.70.+s,25.43.+t,13.75.Cs,13.85.Fb,13.88.te 



I. INTRODUCTION 

Measurements have been made of the pp — > AA reac- 
tion on a transversely polarized frozen-spin target. An 
analysis is presented of data taken at the Low Energy An- 
tiproton Ring (LEAR) at CERN, Geneva at p p = 1.637 
GeV/c corresponding to a center-of-mass energy which is 
78 MeV above A A threshold. 

This experiment, carried out by the PS185 collabora- 
tion, expands upon a series of YY production and re- 
lated experiments Q, 0, 0, which have previously been 
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performed by this collaboration using the same detec- 
tor system. These covered wide-ranging kinematics from 
very near threshold to higher energies at which the larger 
cross section allowed high-statistics studies to be made. 

Spin observables have long been of interest in 
strangeness-production reactions in part because of the 
characteristic strong polarization produced. The hy- 
perons in the final state (the word 'hyperon' will be 
used loosely to include anti-hyperons) lend themselves 
to study of spin dynamics because of the self-analyzing 
power of mesonic hyperon decay. Thus, even before the 
introduction of initial-state polarization in these studies, 
the PS185 collaboration has been able to greatly expand 
the world's supply of data on spin observables, especially 
for pp — > AA. In particular, final-state polarization com- 
ponents of both hyperons and correlations of all spin com- 
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ponents of A and A were accurately measured in addition 
to differential cross sections. 

The wealth of data produced by PS185 ex- 
cited a great deal of theoretical activity |E1E13)1SM0| 
with several models enjoying reasonable success in fitting 
the observations. Two distinct theoretical approaches, 
meson-exchange [MEX] and quark-gluon [QG] inspired 
models, have been used successfully to fit the pp —> AA 
data despite being based on fundamentally different re- 
action dynamics. 

Several authors 0,SI3 have constructed models based 
on t-channel exchange of strange mesons. In order to 
match the observed data, these meson-exchange models 
require a strong tensor interaction and so a spin-flip from 
the initial spin-triplet pp pair to the final spin-triplet AA 
pair. Initial- and final-state interactions are modeled as 
well but are found not to qualitatively change the spin- 
flip character of these models. 

The alternate quark-gluon inspired models 0, EJ are 
based upon an assumed s-channel interaction between 
qq pairs leading to the transformation to an ss pair. In 
existing models the qq pair is assumed to have vacuum 
quantum numbers (0 + for 3 P qq pairs) or gluon quan- 
tum numbers (1~ for 3 Si qq pairs). Again there is a 
dominance of the spin-triplet state (inspired by the em- 
pirical fact that the singlet fraction is, on average, small). 
In fact this construction assures pure triplet transitions. 
(Since the spin of the A reflects the spin of the s-quark, 
in the constituent quark model, triplet ss pairs guarantee 
a triplet AA final state.) Here the assumed fundamental 
s-channel qq — > ss process does not flip the spin at all but 
some spin-flip strength is introduced with the inclusion 
of initial- and final-state interactions. 

The difference in spin-flip predictions of MEX and QG 
models was suggested 01 as a means of distinguishing 
the validity of the two classes of models. Since each model 
predicted strong triplet interactions, final-state spin cor- 
relations were similar and measurement of final-state spin 
alone was not sufficient to distinguish them. Sensitivity 
could be found in measurement of the correlation between 
initial-state spin and final-state spin. In particular, the 
normal-to-normal depolarization and spin transfer were 
expected to strongly select between the two classes of 
models. These observables, often denoted D nn and K nn , 
respectively, are sensitive to the transfer of spin from the 
initial-state proton to the final-state A (in the case of 
D nn ) or A (in the case of K nn ). The 'n' subscripts indi- 
cate that the component of spin considered in each case 
is the component normal to the scattering plane. Since a 
large number of spin correlations will be discussed here, a 
more general notation will be introduced below. In par- 
ticular, D nn and K nn will be denoted as Q[n p ^nj\\ and 
<3[n p , ro^h respectively. 

Experimental study of these particularly interesting 
observables required a polarized target. A frozen-spin 
target was constructed with such small dimensions that 
AA pairs could exit the target before decaying. The suc- 
cess of the present measurements relied upon the PS185 



apparatus, this new target, and the superb properties of 
the LEAR beam. 

It was also noted that such data contained enough 
information to permit determination of not just a few 
observables but the entire spin-structure of the reaction. 
From this, all possible spin observables could then be 
determined. That analysis has been successfully carried 
out and is reported here. A previous publication 
[l4| already reported the measured spin-transfer and de- 
polarizations, extracted using the techniques which are 
explained in the present publication. Interestingly, the 
measured values disagree strongly with the predictions 
from both classes of models, leaving opportunity for fur- 
ther theoretical study. These results are also included in 
the present paper. 



II. SPIN CORRELATIONS 

The density-matrix formalism lends itself to analysis of 
systems composed of ensembles of non-interfering states, 
such as occur in polarized systems. This formalism will 
be used therefore to precisely define the spin-correlations 
and to relate them to the observed final-state distribution 
on the one hand and to the spin-scattering matrix, which 
parameterizes the transition, on the other. 

In this formalism, the expectation value of an ob- 
servable represented by operator O can be written as 
< O >= Tr(Op AA ) where p AA is a 4 x 4 density matrix 
representing the final-state spin information. In Sec. Ill Al 
we write p AA in terms of the initial-state polarization 
vectors, P p and P p , and the pp — > AA transition opera- 
tor M(Q c .m.)- Since the experiment actually detects the 
final-state protons, antiprotons, and pions from the hy- 
peron decays, the formalism is extended in Sec. Ill Bl to 
relate the directions of these particles' momenta to the 
initial-state polarization vectors and the complete set of 
the reaction's spin-observables. 



A. Spin Dynamics of pp — > AA 

The initial-state density matrix for spin-i particles 
with polarization P can be written as 

P =\(l + P.a) 

where X is the identity matrix and Oi are the Pauli ma- 
trices. In the following discussion, notation is greatly 
simplified by defining Pq — 1 and do = 1. Then 

k=Q 
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For an initial state of a proton and anti-proton the initial- 
state density matrix is 

P* = ^ = jE pf°fH°k (i) 

j,k=0 

where <t^ and <7j^ operate in the separate spin-space of 
the anti-proton and proton, respectively, and a direct 
product is implied. 

Meanwhile, the density matrix after the interaction, 
p AA , can be written as an arbitrary linear combination 
of direct products of a B and a B since these span the 
space of Hermitian 4x4 operators 

P KA = £ Wr4o* (2) 
i±,i> — o 

where the notation of a B and a B from Eq. (|TJ have been 
used in anticipation of the fact that an identification will 
be made between the p and A anti-baryon spin spaces 
and between the p and A baryon spin spaces. 

Spin correlations are meaningful observables as long 
as the coordinate system associated with each particle is 
well defined. It is not necessary that the coordinate axes 
used for one particle be aligned parallel to those used 
for another. The following discussion of spin correlations 
applies for any set of coordinates in which the spin com- 
ponents of each particle are measured in some set of axes 
defined in the rest frame of that particle. The specific 
choice of coordinate axes for this analysis is described in 
Sec. Urn below. 

It then follows (since Tria^ay) — 25^) that 

4WV = Tr(a B a B p AA ). 
Substituting this into Eq. gives 
1 3 - 

AA 1 \ " rp / B B AA\ B B fo\ 

p = 4 2^ 1 r ^ a » p > a v °v ■ y s ) 

On the right hand side p AA can be rewritten in terms 
of pP p , which is given in Eq. (JJJ. If M is the transition 
operator from the pp to AA state, 

1 3 

p AA = MpPPA'P = PjPkMafagM^. (4) 

i,fc=o 

Substituting this in Eq. gives the identity 

iU,i/=0 j,k=0 

(5) 

Only the first term contributes to the unpolarized dif- 
ferential cross section, 

= Tr (p AA ) = 1 r r MM U 
Tr(pPP) 4 V ; 



Factoring this out of the sum in Eq. J5) and defining 

Q\jp, k p , p A , u A ] = -^-Tr(a B a B Ma B c7 B M^ (6) 
allows Eq. (J5J to be written as 

3 3 

/^^E J2PfP p k Q[j P ,k p ,p A ^ A }a B a B . (7) 

fi,u=0 j,k=0 

The spin dynamics of the pp — ► AA reaction at any 
production angle, O c . m ., is entirely contained within the 
quantities Q[j p , k p , (X A , v A \. These Q's are the observables 
of the present experiment, which we shall refer to as the 
'spin correlations'. The unpolarized differential cross sec- 
tion is parameterized by Iq. Although the functional de- 
pendence will not be written explicitly, it is to be under- 
stood that Iq and the Q's are functions of O c . m .. The Q's 
would be directly measurable if final-state spins could be 
measured directly (and initial-state polarizations could 
be chosen arbitrarily). For example, measurement of the 
mean value of the product of the ^-component of A spin 
and the ^-component of A spin (given the initial p po- 
larized in the r direction and proton polarized in the s 
direction) would be directly related to Q[rp, s p , r/ A , £a] by 

§ B Tr(p AA a B a B ) _ PfPfQ[r p ,s p ,r, K ,U] 
< ^ £ >- Tr{pAA) ~ Jf 

where the normalization factor is given by 

N = l + PPQ{r p ,0,0,0} + PPQ[Q,s p ,Q,0) 
+ P?P£Q[r p ,s p ,0,0] 

The redundant particle-identification subscripts have 
been introduced in the index list of Q to allow suppres- 
sion of vanishing elements. Subsequently terms such as 
Q[0, s p , 0, 0] will be written simply as Q[s p ]. 

B. Angular Distribution of Decay Products 

In fact AA final-state spin information is not directly 
measured in this experiment. It can be inferred, how- 
ever, from the angular distribution of the decay products 
because of the self-analyzing nature of the A and A de- 
cays. For a A with polarization vector P A , the angular 
distribution of the decay proton in the A rest frame is 
given by 

I p {kP) = l-(l + aP A -k p ) 

where a = 0.642 ± .013 is the self-analyzing power of 
A — > p7r~ and k p is a unit vector in the direction of the 
proton's momentum (in the A rest frame). Similarly, for 

A — * pK + , 

I p (k p ) — -j— (1 + aP A ■ kP) 
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where a = —a by CP-conservation. 

The transition operator T A , representing A — > pir~ de- 
cay must give this observed angular distribution from 

I p (k p ) = Tr(T AP A Tl) = Tr(T A ±(I + P A ■ <?)T A ) 

It then follows that 

Tr(T A a t T{) = ^ if i e{l, 2, 3} (8) 

Z7T 

(where kf is the i'th directional cosine of the proton's 



IfinaliQcm., $c.m.,k p , k p ) = Tr(T A T A p AA Tl T A ) = 



where fc P h as similarly been extended by defining h$ = i . 

Eq. (|10fl provides the connection between the ob- 
served distribution, I final-, and the spin correlations, 
Q\jp, k p , /J, A , v A ], which are functions of the production 
angle c . m . and of the beam energy. If no further con- 
straints were known, the spin correlations could be ex- 
tracted (for each bin of c . m .) from the observed de- 
cay angular distributions. In fact, as discussed in the 
next section, the structure of the spin-scattering matrix 
is subject to constraints which enforce parity and charge 
conjugation symmetry. As a result the spin correlations 
are not independent functions. Rather than attempting 
to determine spin correlations independently, it is prefer- 
able to determine the parameters of the spin-scattering 
matrix and to use them to determine the spin correla- 
tions of interest. This results in improved precision of the 
resulting spin correlations and guarantees that they are 
mutually consistent, obeying all constraints implied by 
the structure of the allowed spin-scattering matrix |l6j . 
It has been shown in reference [3 that measurements 
of pp — > AA — > pTr + pTT~ with an unpolarized beam on 
a transversely polarized target provide enough informa- 
tion to fully constrain the spin-scattering matrix. This 
then allows determination of all spin correlations, even 
those with non- vanishing jp, for example, whose direct 
measurement would require the use of polarized beam. 

C. Coordinate systems 

The emphasis of this experiment was the determina- 
tion of correlations betweens spins, particularly between 
initial-state target spin and and final-state spins. Com- 
parison of initial- and final-state spin is complicated by 
the fact that initial-state particles are moving at relativis- 
tic velocity (/3=.58) in the center of mass. Comparing all 
spins in a common reference frame would then require 



momentum) while 

Tr{T K(H Tl) = Tr(T A Tl) = ±- if i=0 (9) 

Z7T 

For notational simplification, Eqs. (JSJ and @ can be 
combined into a single equation of the form of Eq. © 
by extending the definition of k p , defining feg = — . Then 
Eq. (JBJ) can be used to find the angular distribution of 
the final-state proton and anti-proton resulting from the 
decay of the A and A, 



-^gFE E^fe,^,MA^A]p;p^ (io) 

u,v=0 j,k=0 

I 

a relativistic boost of the spins. Such a boost of 2-D 
spinors is not well defined. Consistent transformation to 
a common reference frame would require use of relativis- 
tic 4-spinors. This complication is avoided by defining 
all such correlations in terms of the spin projections of 
each particle in its own rest frame. As mentioned above, 
it is not necessary for the coordinate axes used to de- 
scribe the spin of one particle be parallel to those used 
to describe the spin of another. Indeed, since there is no 
common boost direction between particle rest frames it 
is not generally possible to define the coordinate axes to 
be all mutually parallel. It is conventional to use helicity- 
based coordinates in which one axis is aligned with the 
particle's center-of-mass momentum direction. 

The final-state polarization direction is inferred from 
the angular distributions of the decay products in the rest 
frames of the A and A. It is therefore natural to express 
spin information for each of these hyperons with respect 
to coordinate axes which are defined in their respective 
rest frames. Similarly, it is natural to express the target 
spin information with respect to axes which are at rest in 
the lab frame. For notational completeness, although the 
beam is unpolarized, its spin is expressed in a coordinate 
system moving with the incident p. 

The specific coordinate systems used in this analysis 
are represented in Fig. ^ All four coordinate systems 
share a common n direction, which is the normal to the 
scattering plane defined as a unit vector in the direction 
of p p x pa (or equivalently in the direction of pp x p^). 
For the coordinate system associated with each particle 
a second axis (l p , lp, 1a, 1a) is defined along the di- 
rection of the particle's momentum with respect to the 
center of mass of the pp system. Finally the third axis 
(rhp, rhp, iSia, hia) is defined for each coordinate sys- 
tem as rhj = ri; x 1^ so each (1,, m^, hi) defines a right- 
handed orthonormal set of basis vectors in the rest frame 
of particle i. 
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FIG. 1: The coordinate axes used to decompose each of the 
individual particles' spin directions are shown. The h direc- 
tion is common. The target polarization direction, which is 
perpendicular to l p is also shown for a typical event in which 
the normal to the scattering plane is at angle $ relative to 
the polarization direction. 

It is of course important to take into account these 
coordinate definitions when the results presented here 
are compared to predictions or to other measurements. 
The meaning of a specific spin correlation depends upon 
the axes with respect to which it is defined. For exam- 
ple, with this choice of coordinates, the angle between 
rhp and riiA depends upon the scattering angle. This 
must not be neglected when interpreting a result such 
as the correlation between the m p -component of the ini- 



tial proton spin and the riiA-component of the final A 
spin. These coordinates were also used in extracting the 
previously published results 0] from this measurement. 
Since those results involved correlations between normal 
components of spin, there was no potential for ambiguity 
in interpretation of those results. 

With these coordinate definitions, n is an axial vector 
while m and 1 are polar vectors for each of the coordinate 
systems. Then for any of the spins expressed in its respec- 
tive coordinate system, a n is scalar under parity inversion 
while oi and a m are pseudoscalar. Parity conservation 
then requires that Q[rp, s pi rj A , £a] must vanish if an odd 
number of l's and m's appear in {r,s,77,£}. Since the tar- 
get is transversely polarized in this experiment Pf is zero. 
This reduces the sum over k in Eq. (|10|l to run only over 
k e {0,m,n}. Furthermore, the unpolarized beam means 
that only the j=0 terms survive. Further constraints on 
Eq. 111)11 result from C-parity conservation which requires 
that Q[jp, k p , [i Al ^a] = Q[kp, j p , v A , MaI- Additional sim- 
plification results from Bohr's rule |l7j| which requires, 
for example that Q[n p , m Al m A ] = — Q[n v , l A ,l A ]. Incor- 
porating all these simplifications, Eq. (|1(J|) reduces from 
a sum over 256 terms to 



1 final (®c.m. i ^c.m. : kP , k^j 

/ (e c . m .)/i6^ 2 



i 

Q[n\] 

Q[n A ,n A ] 

Q[m A ,m A ] 

QHM 

Q[m A ,l A ] 

Q[n p ] 

Q[n p ,n A ] 

Q[n p ,n A ] 

Q[m p ,m A ] 

Q[m p ,l A ] 

Q[m p ,m A ] 

Q[m p ,l A ] 

Q[n p ,m A ,m A ] 

Q[n p , m A , I A ] 

Q[n p J A ,m A ] 

Q[m pi m A ,n A ] 

Q[m p ,l A ,n A ] 

Q[m p ,n A ,m A } 

Q[m pi n Al l A ] 



(akP + akP) 



dak? fcP 
aakfkf 



kfk^ 



{P T cos ($ c . m .) + daP T kPkP cos ($ c . m .)) 
aP T k? n cos ($ c . m .) 
aP T ££cos ($ c . m .) 
aP T kP n sin ($ c . m .) 
aP T kfsm ($ c . m .) 
af T ^ sin ($ c . m .) 
aP T k\ sin ($ c . m .) 
aaP T (kP n k p n cos($ c 



aaP T k P n k p cos ($ c . m J 
aaP T kfk p n cos ($ c . m .) 
aaP T kP n kP sin ($ c . m .) 
aaP T kfkP sin($ c . m .) 
aaP T kP n kP m sin ($ c . m .) 
aaP T kPk P sin ($ c . m .) 



- feffcf cos($ c . m .)) 



Eq. contains just 19 spin-correlations along with 

Iq. Since each term manifests a unique angular depen- 
dence these spin-correlations are, in principle, directly 
measurable by fitting Eq. (| 1 1 11 to the angular distribu- 
tion observed in scattering from a transversely polar- 
ized target. Although this is not the technique employed 
here for determination of spin-correlations, these 19 spin- 



correlations will be referred to as the being " directly mea- 
surable" . 
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III. SPIN-SCATTERING MATRIX 

The transition operator M, introduced in Eq. 10}, 
transforms from the space spanned by direct products 
of proton and anti-proton spinors to one spanned by the 
direct products of A and A spinors. Since it includes the 
spin part of the transition, it is called the spin-scattering 
matrix. It can be represented by a 4 x 4 complex ma- 
trix. It is convenient to construct this operator from a 
combination of direct products of 'baryon operators' and 
'anti-baryon operators'. Baryon operators will be defined 
to be those which transform from the proton spin-space 
to the A spin-space while leaving anti-baryon spinors un- 
affected. Conversely, anti-baryon operators will trans- 
form from p spinors into A spinors. This identification 
of proton spinors with A spinors is just a convenience, 
there is no loss of generality since the 16 direct products 
considered will span the entire space of 4 x 4 Hermitian 
matrices. 

Constructing terms with definite symmetry properties 
is simplified by choosing {X B ,af , cr^, a B } as the baryon 
operators and {X B , af , , <7 B } as the anti-baryon oper- 
ators. Although these operators have the same matrix 
representations as the spin operators, they are actually 
transition matrices having the space of proton spinors 
as their domain and the A spinor space as their range. 
For f] e {0, m, n}, each a B transforms eigenstates of 
the ^-component of proton spin to the same eigenstates 
of the ^-component of A spin, multiplied by the eigen- 
value. Similarly, a B maps anti-proton ?y-eigenstates to 
anti-lambda 77-eigenstates. Under parity inversion, the 
behavior of a B differs from that of af and because n 
is an axial vector while 1 and rh are polar vectors. Given 
that they act upon components of spinors and produce 
components of spinors, the <J B, s have the same parity 
properties as the corresponding components of spin, i.e. 
af is scalar under parity inversion while af and af are 
pseudo-scalar. 

Given these properties, the baryon and anti-baryon op- 
erators may be used to construct a complete set of oper- 
ators having good parity and C-parity and spanning the 
direct-product space. A set of such operators is listed in 
Table ^ along with their parity and charge conjugation 
eigenvalues. 

Conservation of parity and C-parity in the pp — > AA re- 
action requires IS] that only the terms with positive par- 
ity and C-parity contribute to the spin-scattering matrix. 
An arbitrary linear combination of the allowed terms can 
be constructed by weighting each term by the coefficients 
give n in Table [I] A conventional parameterization for M 
|l3 is 

i { (a + b)l B l B + (a - b)afaf + (c + d)a B a B 
M= + {c-d)afaf + e{af l l l3 +l B af) 

Since overall phase is unimportant, the six complex pa- 



Operator 


P 


c 


coefficient 


1 B 1 B 


+ 


+ 


a + b 


„B B 


+ 


+ 


a — b 


„B B 


+ 


+ 


c + d 


B B 


+ 


+ 


c — d 


„BrrB | -rB B 


+ 


+ 


e 


„B-rB tB B 


+ 


— 





„B-jB 1 tB „B 


- 


+ 





B-jB T B B 









„BtB 1 tB „B 




+ 





„BtB -tB „b 









„B B 1 „BB 


+ 


+ 


a 


„B B „B „B 


+ 







B B , B B 




+ 





B B B B 
O m O m " n 









„B B 1 „BB 




+ 





„B B „B B 










TABLE I: A complete set of operators spanning the space 
of operators which transform from pp spinors to AA spinors. 
These operators are constructed to have definite parity P and 
charge-conjugation parity C. The P and C values are also 
listed. The final column gives the coefficient by which each 
term is weighted in forming M. 



rameters {a,b,c,d,e,g} can be represented by just 11 real 
parameters for each O c . m .. Specifically, parameter 'a' will 
be chosen to be real and non-negative while the other five 
parameters have real and imaginary parts. These param- 
eters may be determined by an unbinned 11-parameter 
simultaneous fit to the observed production and decay 
angles of reconstructed pp — > AA — > pit + pTr~ events, as 
described in Sec. IV JJl below. Once the parameters of the 
spin-scattering matrix have been determined, a wealth 
of spin-correlations can be calculated by substituting the 
form of M given by Eq. I|I2|I into the definition of the spin 
correlation, Eq. JfJJl. For example, Table ITT1 gives the re- 
sults, in terms of {a,b,c,d,e,g}, for the 24 spin correlations 
which could be directly measured with a transversely po- 
larized target and unpolarized beam. In fact, Eq. © can 
also be used to calculate other spin correlations which 
would be directly measurable only with longitudinal tar- 
get polarization and/or polarized beam. 

Table m also lists names which have traditionally been 
used to identify some of these spin correlations, A for 
scattering asymmetry, P for final-state polarization, D 
for depolarization, K for spin transfer, and C for final- 
state spin correlations. Here C is also used for 3-spin 
correlations between the target and final state. Caution 
should be used in identifying these spin correlations with 
those in other publications having the same traditional 
name. The precise meaning of spin correlations involving 
m and 1 depend critically upon the choice of coordinates. 
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Spin Correlation 



Traditional 



Jo X Q ( = \Tr(a^a^Mafa^M^ ) 

r _ 1 ri„|2 i l;,|2 , |„|2 , |j|2 , |i2 , |„|: 



Q[0 P ,0 P ,0 A ,0 A ] = 1 

Q[n A ,n A ] C„„ 

Q[n p ,n A ] Ann 

Q[n p ,n A ] 

Q[n A ] = Q[n A ] P n = P„ 

Q[n p ] = Q[n p ,n A ,n A ] 4„ = C nnn 

Q[m A , I a] = Q[Z A , m A ] C mi = C im 
Q[n p ,m A ,m A ] = — Q[n p ,Z A ,Z A ] C„ mm = — C n i, 
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TABLE II: The first column lists the spin correlations which are directly accessible with a transversely polarized target and 
unpolarized beam. The third column gives the indicated correlation (multiplied by 7o) in terms of the parameters of the 
spin-scattering matrix given in Eq. il l'it . The second column lists alternate names for the spin correlations. 



IV. APPARATUS 

With the exception of the polarized target and asso- 
ciated trigger detectors, most of the equipment was the 
same as that used in previous versions of the PS185 ex- 
periment P, 0, 0, 3 ■ A- schematic view of the appara- 
tus is shown in Fig. [21 A compact set of tracking de- 
tectors was located just downstream of the target area. 
Because of the forward boost of the AA system, the tra- 
jectories of both particles passed through these tracking 
chambers. This resulted in a large acceptance for the 
full reconstruction of the charged tracks resulting from 
pp — > AA — » pTt + pn~. Accurate measurement of the 
topology of these charged tracks is sufficient, apart from 
ambiguity of A vs. A, to completely reconstruct the kine- 
matics of each event including the hyperon decay angles 
whose distributions yield information on the final-state 
spin. The ambiguity of A vs. A is resolved by a solenoid 
magnet further downstream which bends the trajectories 
enough to allow the sign of the charge of tracks to be 
determined. 

The tracking detectors consisted of two detector stacks. 
The first was 10 planes of multi-wire proportional cham- 
bers with planes alternately oriented at ±45° relative 
to the horizontal (the U and V directions). These had 
a pitch of 1.27 mm and a separation of 1 cm between 
planes. The second stack was 13 planes of drift chambers 
oriented vertically and horizontally (the X and Y direc- 



tions) with an average separation between planes of 1.35 
cm. The 4 cm drift cells were read out by a pair of sense 
wires separated by .42 mm, which resolved the usual left- 
right ambiguity. A set of three similar but larger drift 
chambers was located inside the solenoid magnet to de- 
termine the direction of deflection of tracks by the 1 kG 
magnetic field. 

Since the beam passed through all these chambers, it 
was necessary to desensitize the center of each chamber. 
This was done by electro-plating additional metal onto a 
3 to 10 mm length of the relevant sense wires to thicken 
them at the position at which the beam would pass, pre- 
venting gas amplification. 

Four planes of microstrips with 100 fim pitch were lo- 
cated upstream of the target. These provided tracking of 
the individual incident p for a fraction of the events. Four 
planes provided no redundancy so a beam track could be 
reconstructed only for events in which exactly one clus- 
ter was found on each plane. This applied for only about 
45% of events because of beam pile- up and detector-aging 
in the high intensity beam. For other events the average 
beam track was used. 

The target was enclosed within a 4.2 cm diameter vac- 
uum vessel. Scintillation counters were used to form a 
trigger which exploited the charged-neutral-charged sig- 
nature of pp — * AA — > pir + pir~ events by selecting events 
in which a p entered the target vessel, no charged parti- 
cles left the target region and at least one charged particle 
exited the tracking chambers. Figure ^jp shows the scin- 
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FIG. 2: a) A schematic plan view of the detectors, magnet, and target. The solenoid field, shown in the positive y-direction, 
could be reversed to reduce systematic uncertainties. Tracks from a typical event are superimposed, b) An expanded view of 
the target region shows the beam-defining scintillators and the veto scintillators used to select events in which only neutral 
particles exit the target. The hodoscope planes completed the trigger by requiring charged particles downstream of the decay 
region. 



tillators used to require an incident p and to veto events 
in which charged particles exit the target area. Figure 
|2t shows the two scintillator hodoscope planes which in- 
dicated a charged decay product. 

The most significant change in the apparatus from pre- 
vious versions of PS185 was the change from a polyethy- 
lene active target to a transversely polarized frozen spin 
target|2(|. As described above, it has been shown [T3I 
that this greatly expands the accessible spin information. 
The target itself was a solid 9 mm long, 6 mm diam- 
eter cylinder of frozen Butanol (C4H9OH) doped with 
TEMPO (C 9 Hi 8 NO)|23. The axis of the cylinder lay 
along the beam direction. Figure 01 shows the nested 
cryostats and the windows which allowed the beam to 
enter and the AA to exit. An extremely compact de- 
sign was achieved to hold the target at 60 mK within 
a room-temperature outer vacuum vessel of only 2.1 cm 
radius. Such a compact design was critical to the success 
of the experiment as it allowed a significant fraction of 
the pp — > AA — > pw + pir~ events of interest to be selected 
by their neutral intermediate state. With a large target 
vessel the hyperons would have decayed internally and 
the trigger information would have been lost. 

The target was polarized roughly every 22 hours by 
surrounding it by a superconducting solenoid, with a field 
of up to 5 T, and pumping with microwaves to achieve 
dynamic nuclear polarization. Initial polarizations up 
to 75.3% were achieved. In frozen spin mode, with a 
holding field supplied by an internal solenoid, polariza- 
tion lifetime was roughly 100 hours with beam on target. 
The average polarization over the run was about 62%. 



The polarization direction was vertical (transverse to the 
beam) and could be chosen to be positive or negative to 
reduce systematic errors. 

The beam rate on target was approimately 6 x 10 5 an- 
tiprotons per second. The integrated beam on target for 
this measurment was 1.6 x 10 11 antiprotons. Scaling by 
the data acquisition live time of 79.2% yields an effective 
integrated intensity of 1.27 x 10 11 antiprotons. 

V. DATA ANALYSIS AND CROSS SECTION 
RESULTS 

A. Event Reconstruction and Kinematic Fitting 

The events of interest, pp — > A A — > pTt + pTr~ , are char- 
acterized by their two-'Vee' structure resulting from the 
decay of the neutral hyperons. The first goal of the 
analysis is to extract the small number of candidates 
for this event topology from the very large number of 
events (largely hn followed by n annihilation) which sat- 
isfy the charged-neutral-charged trigger condition. Less 
than 0.1% of the recorded events were ultimately found 
to be consistent with the pp — > AA — > p~K^pn~ hypothe- 
sis. The next goal is to determine which candidate events 
match the constraints enforced by energy and momentum 
conservation on real pp — > AA — > pTr + pTr~ events. This 
is done by kinematic fitting of the event topologies in 
the tracking chambers, which are in a region free of mag- 
netic field. This fit also provides the best estimate of the 
production and decay angles of interest and of their cor- 
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FIG. 3: (Color online) Schematic cross-sectional view 
through the cryostat system (four concentric vertical cylin- 
ders) holding the frozen spin target (a horizontal cylinder). 
Coordinates are indicated relative to the center of the target. 
The antiproton beam entered through the four thin windows 
on the left and interacted in the frozen butanol target. Hy- 
perons escaped with relatively little interaction through the 
large exit windows shown on the right. Dashed lines represent 
20 fj,m titanium windows while dash-dot lines represent 40 /im 
aluminum windows. Coils which provide the holding field are 
shown shaded. The trigger selected only those events in which 
the hyperons survived long enough to escape the cryostat and 
pass the veto detectors before decaying. 



relations. Additionally, information from the three drift 
chambers within the solenoidal magnetic field is used to 
determine the sign of the charged particles. If at least one 
particle's sign can be clearly determined, then the A-A 
ambiguity can be resolved and the event can be used in 
determining spin correlations. 

When possible, microstrip information was used in the 
kinematic fit to define the position and direction of the 
incoming p. When this was not possible, average beam 
position and direction were used. Because of the very 
small emittance of LEAR's adiabatically cooled beam, 
this was sufficiently precise. 

Drift chamber wire positions and time-to-distance cali- 
brations were determined empirically from straight track 
data. This allowed drift chamber hit positions to be ac- 
curately determined, typically with better than 200 /im 
resolution. 

For candidate pp — > AA — > p7T + p7r _ events, hit posi- 
tions were determined from MWPC wire positions and 
from drift chamber wire position and time. A significant 
walk correction was made for drift chamber hits based 
upon time-above-threshold. Hits from all planes sharing 



a common readout direction (X,Y,U, or V) were then 
searched for 2-D track projections (X-Z, Y-Z, U-Z, or 
V-Z). Drift chamber hit positions on a 2-D track candi- 
date were iteratively improved by correcting drift-time 
to position conversion to reflect the apparent slope of 
the track. To allow for crossing tracks, a hit could be 
included in more than one 2-D track. A maximum track 
angle of 60° was allowed in each projection since tracking 
was inaccurate beyond that range. Losses due to that cut 
were accounted for as acceptance losses as discussed in 
the Scc.lVDl 

Combinations of three or more of these 2-D tracks were 
then considered as candidate 3-D tracks. The 45° rota- 
tion of the MWPC projections relative to the drift cham- 
ber projections helped eliminate spurious combinations. 
Confidence-level cuts based on \ 2 were used to determine 
whether sets of 2-D projections were consistent with a 
3-D track. Tracks constructed from just two 2-D projec- 
tions were considered only if neither projection could be 
used to form a 3-D track having more projections. 

Candidate Vee's were formed from pairs of 3-D tracks 
having a distance-of-closest-approach consistent with 
zero. These Vee's were rejected if they were not consis- 
tent with in-flight decay of a A in the momentum range 
(471-1161 MeV/c) expected for pp — > AA — > pir + pw~. 
A single 3-D track was allowed to be included in more 
than one Vee to avoid losing a real Vee by having a track 
mis-assigned to a false Vee. 

Pairs of Vee candidates (which did not share any 3-D 
tracks) were then considered as candidates for kinematic 
fitting. The pair was first tested for rough consistency 
with the kinematic hypothesis, (e.g. Transverse com- 
ponents of momentum, calculated independently from 
the topology of each Vee, should be opposite and equal, 
within errors.) 

For all tracking cuts, simulated events (discussed in 
Sec. IV B|) were used as a guide in setting confidence- 
level cuts to avoid cutting good events. The simula- 
tion included effects due to finite resolution and multiple- 
scattering. 

Kinematic fitting was based upon the fact that the 
ideal topology of an event (neglecting finite resolution, 
multiple-scattering, and interactions) can be completely 
predicted in terms of 14 parameters: 

— 3 components of beam momentum 

— 3 coordinates of production vertex 

— 2 decay lengths 

— 2 production angles (0 c . m . and <fr c .m.) 
4 decay angles (0a, </>a, a , 4>a )■ 

The production and decay angles are of greatest interest 
since they hold the information on the spin dynamics. By 
varying these 14 parameters it should be possible to find 
a hypothesis which is consistent with the hits on the ob- 
served Vee-pair, within errors. If satisfactory consistency 
cannot be achieved then the pair can be rejected as not 
originating from pp — » AA — > pir + pir~ . If an acceptable 
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fit can be found, then it gives the best estimate of the 14 
parameters. 

When evaluating consistency of a hypothesized set of 
parameters not only the measured track hits were used 
but also a data point was included to represent, with ap- 
propriate errors, the knowledge of the beam energy and 
the beam direction for that event. Individual hits (as op- 
posed to tracks) were treated as measurements in evalu- 
ating the goodness of fit. The errors on the measurements 
were not treated as being independent, however. Evalu- 
ation of goodness of fit took into account the fact that 
the errors in hit positions on any given track were cor- 
related because of multiple-scattering. This correlation 
was increased by the fact that the scattering did not al- 
ways happen uniformly along a track but could be greatly 
increased at points where the particle hit a wire of the 
tracking chambers. Simulated events were used to study 
the covariance introduced by multiple-scattering and to 
ensure that the fitting procedure properly accounted for 
it. 

A generalization 0] of the Levenberg-Marquardt 23J 
method was used to adjust the 14 parameters to minimize 
a likelihood statistic (analogous to x 2 ) which accounted 
for covariance due to multiple-scattering. The estimated 
errors on parameters were assigned based on the covari- 
ance matrix from the fit. As an example of the accuracy 
of event-reconstruction, the azimuthal production angle 
$ c . m was typically determined with an r.m.s. error of 
less than 0.4° (except when it diverged near the poles at 
sin(0 c . m .) = 0). Also cos(© c . m .) had a mean r.m.s. error 
of roughly 0.04 for the worst-case events (0 c . m . w 90°) 
falling roughly linearly to less than 0.004 for 0° and 180° 
scattering. 

The fit was heavily over-constrained by the measured 
topology of the event, along with the 12 constraints due 
to 4-momcntum conservation at all three vertices. A 
rough measure of fit-quality, Q, was calculated by treat- 
ing the likelihood statistic as if it were ^-distributed and 
calculating the confidence level. A flat distribution over 
< Q < 1 would be expected from a true \ 2 statistic. 
The actual fit-quality distribution has a large peak be- 
low Q — .008 resulting from unrelated background events 
(such as anti-neutron annihilation) which produce large 
numbers of tracks. This very sharp peak was cut early 
in the analysis and will not be included in the discussion 
of fit-quality which follows. The remaining Q distribu- 
tion is shown in Figure 0] The peak at high fit-quality 
(0.9 < Q < 1.0) is understood as resulting from events 
in which no track was substantially deflected by hitting 
a chamber wire. Conversely, the peak a low Q results 
largely from events in which one or more track suffered 
significant deflection. As shown in Figure^] Monte Carlo 
simulation correctly predicted the general shape of the 
observed distribution, including the skew towards high Q 
and the peaks at large and small Q. The small residual 
differences are believed to result from imperfect modeling 
of the multiple scattering and matter distribution. 

Events with Q < 0.1 were rejected. This cut not only 
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FIG. 4: Distribution of fit-quality, Q, found by treating lik- 
lihood as if it were distributed with a \ 2 statistic. Shown 
dashed is Monte Carlo prediction, normalized to the data over 
the interval .1 < Q < 1. 

cleanly rejected the spike of background events at very 
low Q, it also rejected a significant fraction of the events 
from quasi-free AA production on a bound proton in a 
carbon, oxygen, or nitrogen nucleus in the target. The 
efficiency of rejection of such events was tested by fitting 
the pp — > AA — > pir + pTr~ hypothesis to a dataset which 
was collected using a pure carbon target. These events 
showed a peak at low Q, with about half the events being 
eliminated when events with Q < 0.1 were cut. The re- 
maining tail extended across the Q distribution. Scaling 
that data to the number of non-Hydrogen nuclei in the 
doped Butanol target, the quasifrcc contamination can 
be estimated as 1.7±0.1% of the final AA data set. Since 
the protons involved are unpolarized, this is expected to 
appear mostly as a dilution of the extracted spin observ- 
ables. This effect has been included in estimates of the 
systematic errors. 

The cut was expected to reject approximately 9% of 
good pp — > AA — > pTT + pir~ events and so would be ex- 
pected to not significantly bias spin correlations unless 
there was a strong angular variation in the losses. Monte 
Carlo simulated distributions for azimuthal production 
angle and individual decay angles matched those of the 
actual accepted angular distributions, without indication 
of angular variation of the losses. Furthermore, analysis 
of Monte Carlo simulated data with the same analysis 
cuts showed the reconstructed spin-correlation and spin 
scattering matrix parameters to be in agreement, within 
expected statistical errors, with those used to generate 
the Monte Carlo events. 

Kinematic fitting in the field-free region cannot distin- 
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guish the A from the A. Tracks were extended into the 
solenoid and the three drift chambers there were searched 
for triplets of hits which could be used to assign a charge 
to at least one of the tracks. The hits were compared 
to predictions based upon the reconstructed momentum 
and angle of the track and both possible charges. Ex- 
pected covariance of the hits, which was large because 
of multiple-scattering in the coils of the solenoid, was 
taken into account. Any track which was well fit by one 
charge hypothesis acted as a 'vote' for which Vee was 
the A. In principle a single vote was sufficient to re- 
solve the ambiguity and properly identify all four tracks. 
Events with multiple votes occasionally had two tracks 
which voted differently on which Vee should be identified 
as the A. Additionally, 8% of the events were unusable 
because they had no vote or had conflicting votes. 

A naive estimate of the error rate of these individual 
votes can be achieved by assuming the error rate is a con- 
stant, uncorrelated with the number of votes. The rate of 
inconsistency between the votes in two- vote events would 
then imply a 3.1% error rate on individual tracks. As- 
suming that the same error rate applies for those events 
which have only a single vote leads to the estimate that 
in total the A — A identification was interchanged in 1.1% 
of the reconstructed events. 

More careful evaluation however shows that this naive 
estimate under-predicts the actual rate of inconsistency 
for multiple- vote events, indicating a correlation between 
number of votes and error rate. Without adjustable pa- 
rameters, the Monte Carlo simulation described in the 
next section makes an excellent prediction of both the 
observed distribution of number of votes per event and 
the rate of inconsistent votes as a function of the number 
of votes. This gives some confidence in the Monte Carlo 
prediction that the actual fraction of events in which the 
A — A identification was interchanged is only 0.7% with 
no marked dependence on O c . m .. The smaller value re- 
sults from the Monte Carlo's prediction of a lower error 
rate for single- vote events than for two- vote events. 

Combining these two estimates, a (0.9±0.2)% misiden- 
tification rate was assumed when correcting for contam- 
ination as explained in Sec. IV Dl The estimated error 
of ±0.2% was included in the systematic error analysis. 
This contamination was negligible in all but the two most 
back-angle bins. 

A total of 30818 events were successfully kinematically 
reconstructed as AA events with the A-A ambiguity re- 
solved. 



B. Monte Carlo Simulation 

An understanding of the angular dependence of accep- 
tance is critical to successful extraction of spin scattering 
information from angular distributions. The acceptance 
function used to extract spin observables was evaluated 
using a simulation designed to incorporate empirical de- 
tector response along with predicted particle interactions. 



The Monte Carlo simulation was also used in tuning algo- 
rithms and setting cuts to optimize tracking, Vee-finding, 
and fitting. Simulated events were also used to study sys- 
tematic errors. 

The GEANT-based simulation included multiple- 
scattering, (5-ray production, and hadronic interactions. 
The latter was especially important for the final-state 
p, which has a large annihilation cross-section. The de- 
scription of the mass distribution included the target, 
cryostat, scintillators, chamber foils, gas, and wires. 

The position and response of detectors were deter- 
mined empirically when possible and used as input for 
simulation. Positions of trigger scintillators were deter- 
mined by tracking through the micro-strip detectors us- 
ing data from dedicated calibration runs taken with a 
thick scatterer upstream of the microstrips. These data 
were also used to determine tracking chamber positions 
and to measure their efficiencies as a function of track 
slope and position. 

Simulation of the wire chambers included observed de- 
crease in chamber efficiencies near the sense wires, the 
observed effective size of the 'dead spot' built into the 
center of each plane, and the observed decrease in effi- 
ciency on the neighboring drift chamber sense wire due 
to field distortion at the 'dead spot'. 

The effect of the trigger scintillators was an important 
component of the simulation. Use of trigger scintillators 
was essential to select out the rare events of interest. The 
cryogenic target however made it impossible to place the 
scintillators as close to the production point as had been 
done in all previous versions of the PS185 experiment. 
Therefore a large fraction of all pp — > AA — > pTr + pir~ 
events were lost because at least one hyperon decay oc- 
curred upstream of the veto scintillators. 

Figure [S] shows a comparison between simulated and 
measured distributions, as an example of a test of the 
simulation. The z component of vertex position of A 
decay is shown. The distribution is seen to be well pre- 
dicted by simulation, including the sharp rise due to the 
position of the veto scintillators. The target, as shown 
in Fig. |2t>, is located upstream of this position but the 
veto scintillators prevent triggers for events in which the 
decay occurs further upstream. 



C. Differential Cross Section 

The $ c . m . -averaged differential cross section for AA 
production at pp — 1.637 GeV/c can be found from the 
spin-matrix parameters {a,b,c,d,e,g}, determined as de- 
scribed in the following sections, as 



da 1 da T 

> = n = h 

dilA 2-7T dcos6 c . m . 

= ±{a 2 + b 2 + c 2 +d 2 + e 2 +g 2 }. 
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FIG. 5: (Color online) Comparison of the hyperon decay ver- 
tex z-coordinate distributions between data (solid line) and 
simulated data set (shaded). 
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FIG. 6: Results for the ^-averaged differential cross section 
(solid) . The previous measurement at 1.642 GeV/c, scaled 
by a factor of 1.26 to match present integrated cross section, 
is superimposed (open diamonds). 



The differential cross section so found is shown as dark 
points in Fig. |SJ The cross section is essentially iden- 
tical to that found by counting AA events within each 
cos(0 cm .) bin and correcting for mean acceptance over 
the angular distributions found below. 

The unequal bin widths seen in Fig. were chosen to 
give roughly equal statistics of observed AA events in 
each of sixteen cos(O c . m .) bins. This allowed a stable fit 
of the eleven spin-matrix parameters to be performed in 
each bin. Uniform binning would have sacrificed forward- 
angle cos(0 c . m .) resolution or resulted in low-statistics 
fits at back angles leading to instability and large er- 
rors. The cross section is discussed here, in advance of 
the explanation of the fitting procedure, to motivate this 
choice of cos(0 c . m .) bins, which applies to all subsequent 
discussion of the matrix parameters. 

The open points in Fig. [fj] show a renormalized ver- 
sion of a previously published result Q from an ear- 
lier version of the PS185 experiment using unpolarizcd 
polyethylene target cells surrounded by scintillator. That 
measurements were made at pp = 1.642 GeV/c which is 



very nearly the same beam momentum as the present 
data, pp = 1.637 GeV/c. For the purpose of compar- 
ison of shapes, the older result has been scaled up by 
a factor of 1.26 to match the integrated cross section 
of 81.1 ± 0.5j^y'g fib measured in the present experiment 
(where the first error is statistical and the second is sys- 
tematic). The systematic errors assigned to the present 
measurement are larger than on earlier ones because the 
cryogenic polarized target introduced larger uncertain- 
ties in target thickness, target position relative to trig- 
ger counters, and hadronic interactions. As described in 
Sec. IV El these systematic errors have been realistically 
estimated so an explanation is required for the apparent 
discrepancy between the present and previous determina- 
tions of total cross section. The normalization discrep- 
ancy in early PS185 results has already been described in 
an earlier publication 3] . Hadronic interactions and mul- 
tiple scattering were not included in the custom-written 
Monte Carlo simulation code used to determine accep- 
tance corrections in early analyses of PS185 data [j], In- 
versions of a GEANT-based simulation have been used 
in calculating acceptances for more recent results 0, 0] • 
Inclusion of hadronic interactions and multiple scatter- 
ing increased the estimated yield by 8-12%. With a 
10 ± 2% adjustment in normalization, the older result 
gives a — 70 .4 ± 0.4 ± 2.2 fib which roughly agrees with 
the present result within errors. The spin-correlations 
are insensitive to any systematic error on overall normal- 
ization. 

Comparison of the two data sets in Fig. shows some 
difference in shape of the distribution at forward angle. 
The present analysis took into account correlated errors 
due to track deflections caused by multiple-scattering. 
This allowed events to be reconstructed which might oth- 
erwise have been lost. The older analysis did not al- 
low for multiple-scattering in event reconstruction. The 
lost events would not have been compensated by accep- 
tance corrections since multiple-scattering was not in- 
cluded in the Monte Carlo simulation. Furthermore, esti- 
mates of expected losses, found by disabling the multiple- 
scattering correlation in fitting in the present analysis, in- 
dicate that the effect is greatest at forward and back an- 
gles, where hyperon momenta are low. The sharper peak 
at forward angles seen in the present data is therefore 
believed to be accurate while the older data is slightly 
distorted by multiple-scattering losses. 



D. Fitting Spin-scattering Matrix Parameters 

For each bin in c .m., Eq. (|llf> represents the dis- 
tribution of events across a 5-dimensional space of an- 



gles ($ c . m ., cos(6» A ), (j) K , cos( 



'A J 



which will be 



represented by a 5-dimensional vector, v for notational 
convenience. With only roughly 2000 events in each 
©cm. bin, performing a simple x 2_n t by subdividing 
each of the five coordinates into bins is excluded. An 
unbinned maximum-likelihood fitting technique was em- 
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ployed. This is often called simply 'maximum-likelihood 
fitting' a name which fails to distinguish from other fit 
methods such as x 2 minimization, which also correspond 
to a maximum likelihood. 

Unbinncd fitting is a limiting case of fitting with Pois- 
son statistics. If the data were binned into K bins with 
n k being the number of events in bin fc, then the likeli- 
hood of the observed data set would be 



K 

C=l[e 

fc=i 



)Ag fcf W^fc.Q) A ^)"* i 

1 n k \ ' 



where 4>{v k > a) is the probability density function which 
depends on the eleven real parameters of the spin- 
scattering matrix, here represented as the 11-dimensional 
vector a. The volume of the k'th 5-dimensional bin is 
Av k . If the number of bins, K, now becomes large com- 
pared to the number of events, N, then n k will be zero 
for most bins and will be unity for TV bins. In this limit 
the likelihood becomes 



K 



N 



£=ne 

k=l 



-tfi(v k ,a)Av k 



]J ( f > {v J ,a)Av j 



where the second product runs over only the occupied 
bins. Neglecting background contamination, <j> is the in- 
tensity distribution, described by Eq. (| 1 1 1) - scaled by the 
integrated luminosity of the experiment, the solid angle 
corresponding to the c m .-bin considered, and the aver- 
age acceptance probability for coordinates v k , v k . 

Taking the log of C introduces a sum of ln(Awj) terms 
which would diverge in the limit of infinitesimal bin size. 
Since these terms are independent of a, discarding them 
does not change the value of a which maximizes the like- 
lihood. Let M! represent — 21n(£) with these terms dis- 
carded. Then in the limit of infinitesimal bin size, 



M' = 2 



dv 



N 



ln(0(U),a)) 



Maximization of the likelihood for a data set {vi 
vn} can be accomplished by finding the pa- 
rameters, a which minimize A4' . Writing <p(v,a) 
as fi(v, a)A(v) where A is the acceptance and fi 
is the remainder of the probability density function 
(which, neglecting background contamination, would be 

Ifinal(Q c .m., $c.m., k p , k p ) scaled by the integrated lumi- 
nosity and the solid angle), 

N 

M' = 2 fi(v,a)A(v) dv - 2 2J In (>(#.,•, «)) 

N 

The last term, which is independent of a, can be dis- 
carded without affecting the position of the minimum. 



This gives the final function to be minimized, which will 
be called M. 

The first term of A4, which incorporates the accep- 
tance function, could in principle be evaluated by Monte- 
Carlo integration (with each simulated event being pro- 
cessed by the analysis routines to see whether it would 
be successfully reconstructed). However, minimization 
of M. would then require prohibitive re-evaluation of this 
term for each new set of parameters a being tested. For- 
tunately the structure of If i na i in Eq. ifTTl) permits a great 
simplification. Each of the twenty terms of Eq. (|llf> can 
be written as a product with an a-dependent dynamic 
term (containing the Q's and Iq) multiplying a purely 
geometric term, Gi(v), which depends on v but not on 
the parameters a. So /i can be written as 

20 

v{v,a) = 22 Di{a)Gi(v) 



allowing the first term of Ai to be simplified since 

/20 „ 
fj,(v,a)A(v) dv = ^A(a) / Gi(v)A(v) dv 

i—l 
20 



where the weights Wi are the moments of the acceptance 
which are independent of a and so need to be evaluated 
only once by Monte-Carlo integration. To perform this 
integration the simulation was used to generate events 
uniformly distributed in w-space without weighting to 
match the observed spin-correlations. The fitting pro- 
cedure is then to search for a which minimizes 
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M = 2 W l D l (a) - 2 HKvj,3)) 



As explained above, an estimated (0.9 ± 0.2)% rate 
of mis-identification of A and A causes contamination 
of & c .m. bins with data from the supplementary angle. 
Because the differential cross section is strongly forward- 
peaked, this is a negligible effect except for the two most 
back-angle bins. The predicted mis-identification rate 
at forward angles leads to an 8% contamination of the 
most back-angle bin. To allow for this, the /i(v) used 
in back angle bins was not simply I final (v, a) scaled 
by the integrated luminosity and solid angle. Rather, 
I final (v, a) was replaced by an appropriately weighted 
linear combination of I final (v, a) and the background 
term Ifinai^reversed: & forward) where v reverse d is the cor- 
responding point in v-space reached by reversing A and A 
identification and df orwart i is the set of parameters which 
applies for the forward-angle bin. Fitting of forward- 
angle bins was carried out first to allow this contamina- 
tion to be correctly modeled when fitting the back-angle 
bins. 
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Minimization in the 11-dimensional a-space was ac- 
complished by the Polak-Ribiere conjugate gradient 
method |2^. To guard against false minima, each min- 
imization was carried out multiple times with different 
randomly-chosen starting points. For most cos(O c . m .) 
bins a common minimum was found in every search. In 
the worst case the fit converged on false minimum less 
than 65% of the time. Additionally, Monte Carlo simu- 
lated data sets were generated with similar statistics and 
spin-correlations to the actual data. Fitting these simu- 
lated data sets demonstrated the robustness of the tech- 
nique for converging on the proper minimum and allowed 
determination of a scale of confidence level for numerical 
values of M. Unlike % 2 , M has no a priori expected 
value for good fits because it is arbitrarily offset from the 
true log-likelihood. There was no indication that M val- 
ues obtained for the real data were systematically higher 
than what was expected for a good fit, based on values 
obtained for corresponding simulated data sets. 



Use of the curvature matrix, based on the assump- 
tion of parabolic behavior of M near the minimum, was 
found to give inaccurate estimates of the errors on the a 
parameters. Errors were instead determined by using a 
'brute-force' search of the space around the minimum to 
find the maximum possible change in each parameter (in 
conjunction with changes in all other parameters) con- 
sistent with an increase of less than one in the value of 
M. Since M differs from — 21n(£) only by a constant 
offset, a given change in M has the same interpretation 
as the equivalent change in — 21n(£). These maximum- 
acceptable (positive- and negative-) changes in each pa- 
rameter will be referred to as the Ttr' errors since they 
cover the confidence interval which would be covered by 
1(7 in the case of Gaussian errors. Similarly, a '2cr' er- 
ror bound on each parameter was found by finding the 
maximum possible (positive or negative) change in each 
parameter for which M would exceed its minimum value 
by less than 4. 



Once the best-fit parameters a have been found, Iq 
and all the Q's can be determined (even those Q's not 
directly measurable in an experiment with unpolarizcd 
beam and a transverse target polarization). But, since 
the errors on the parameters are highly correlated, the 
error on each spin-correlation cannot be determined by 
simple lowest-order error propagation. Even if the corre- 
lation matrix were determined, error propagation would 
be unreliable. A far superior estimate of the error on each 
spin-correlation comes from the same 'brute-force' search 
for AM — 1 and AM = 4 regions. At each point, all 
quantities of interest, such as spin-correlations, were cal- 
culated. The maximum positive and negative excursions 
of each such quantity from its best-fit value, consistent 
with AM < 1 and AM < 4, were taken to be the 'la' 
and '2cr' errors, respectively, on that calculated quantity. 



E. Systematic Errors 

Overall normalization errors would not affect measured 
spin-correlations but would directly change Iq and equiv- 
alently the $ c . m . -averaged differential cross section pre- 
sented in Fig. Since the focus of the experiment was 
on spin-correlations, overall normalization was not con- 
trolled as carefully as other aspects of the experiment. 
The length of the cylindrical frozen Butanol target was 
measured to be 9.0 ±0.5 mm. An upper limit on the mis- 
alignment of the target axis relative to the beam was 20° 
which could increase effective target thickness by up to 
6%. The uncertainty in effective position of veto scintil- 
lators relative to the target was estimated at 1 mm which 
was found, through Monte Carlo simulation, to cause less 
than 4% uncertainty in normalization. Statistical error 
in the estimate of quasi-free contamination introduced a 
negligible systematic error. The overall fractional sys- 
tematic uncertainty, estimated by adding these contribu- 
tions in quadrature, is +7.1%/— 9.3%. This gives the sys- 
tematic error, quoted above, on the total cross section, 
a = 81.1 ± 0.5j]y5/i6. This same fractional systematic 
error applies to the bins of < ^ > in Fig. While 
the normalization error cancels in spin-observables, all 
spin-matrix parameters would scale by the square-root 
of any overall normalization factor. Since this is a com- 
mon factor on all terms, it is not included in the error 
band assigned to each of these parameters. It should 
be remembered that an overall normalization error of 
+3.5%/— 4.8% applies to the spin-matrix parameters. 

The reconstructed angles for each event, Vk — 
{^c.m., cos((9 A ), cos(#a), </>a), are not known ex- 
actly but are extracted with known errors and corre- 
lations. The method of unbinned maximum-likelihood 
fitting treats each event as a precise point in the 5- 
dimensional w-space and does not incorporate a method 
of allowing for the finite errors on these points. Neglect- 
ing these finite errors introduces a source of systematic 
error in addition to the statistical error discussed above. 
The method of estimating these systematic errors was 
based on the use of simulated data with errors and event 
statistics similar to the actual data. These simulated 
data sets were generated based on spin-correlations cho- 
sen to nearly match those found in the data. Simulated 
events were kinematically fit giving Vk 's with errors sim- 
ilar to the real data. Unbinned maximum-likelihood fits 
were then used to extract best-fit values of the spin- 
transfer matrix parameters, which were then used to cal- 
culate spin-correlations and all other variables of interest, 
as described in the next section. The fitting process was 
then repeated using the ideal VkS which had been used 
by the Monte Carlo to generate the simulated data. The 
differences between best-fit values of each variable fit to 
the simulated data and the best-fit values fit to the ideal 
Wit's was a measure of the systematic error in that variable 
because of the finite resolution of the Vk's. Changes in 
each variable of interest were calculated independently, 
rather than using error-propagation to extract the ex- 



15 



pected effect. Ideal fits were compared to fits using kine- 
matically fit finite-resolution data for ten different simu- 
lated data sets. Additionally, since kinematic fitting was 
quite time-consuming, the statistics of this study were 
augmented by generating an additional 20 data sets by 
simply smearing each ideal Vk by its estimated resolution 
rather than simulating an event and operating the full 
analysis chain on it. Estimates of systematic error due to 
finite resolution were thus found as the r.m.s. shift from 
the zero-resolution value for each variable of interest in 
each Q c .m. bin. This contribution typically dominated 
the overall systematic error estimates which are shown 
with each variable in the next section. 

The non-hydrogen nuclei in the target were unpolar- 
ized so quasi-free events mis-identified as free pp — ► 
AA — ► p7r + p7r~ events will exhibit no correlation to target 
polarization direction. The effect of the small quasi-free 
contamination can therefore be estimated by introducing 
an appropriate fraction of isotropically-distributed simu- 
lated events in place of some of the events of a simulated 
data set. The estimated quasi-free contamination (from 
analysis of carbon-target data) was typically only about 
1% but rose to 3% at farthest back-angle. Systematic 
errors for each variable of interest were again found by 
determining the r.m.s. shift of that variable caused by 
the simulated contamination. This small contamination 
generally caused smaller systematic errors than did the 
effects of finite resolution. 

Several precautions were taken to reduce systematic 
effects due to target polarization. The direction of target 
polarization was reversed during data collection to can- 
cel systematic effects due to any up-down asymmetry of 
the detectors. Average target polarization was measured 
by Nuclear Magnetic Resonance measurements at the be- 
ginning and end of each data-collection period (between 
target re-polarizations) to improve estimates of its value 
at intermediate times. Data-collection periods were typ- 
ically limited to a quarter of the polarization lifetime to 
keep the polarization high and to improve accuracy of 
estimated polarization. The probability density function 
used to weight each measurement in the unbinned fit was 
not based upon an average polarization but on the best 
estimate of the polarization at the time that event was 
recorded. 

An additional systematic error results from uncertainty 
in initial polarization and from possible inhomogeneities 
in target relaxation. The fractional error in measured 
target polarization is estimated at AP/P = 2.3%. The 
average polarization measured at the end of a data- 
collection period could not determine whether depolar- 
ization was non-uniform due to beam heating. However 
the target was kept cold enough that relaxation time was 
not strongly temperature-dependent. Inhomogeneity of 
depolarization is estimated to contribute an uncertainty 
of at most AP/P — 3.8% just before repolarization. 
These errors combined to give a worst-case estimate of 
AP/P — 4.5% at the end of the data-collection period. 

The maximum systematic error due to polarization was 



estimated by shifting polarizations by 4.5% and deter- 
mining the size of the shift of extracted observables of 
interest. This was found to be a smaller error than that 
due to resolution or quasifree background. All three ef- 
fects were added in quadrature bin-by-bin for all variables 
of interest to obtain the overall systematic error estimate 
which is plotted as an error band at the bottom of plots 
of each variable of interest. 



VI. SPIN CORRELATION RESULTS 

Figures and |H1 show the best-fit values of the 11 pa- 
rameters of the spin-scattering matrix fit to the data in 
bins of Q c .m.- Parameter 'a' is chosen to be real (and 
non-negative) so Im(a) is zero. On each point the dark 
bar indicates the 'lcr' error range while the light bar rep- 
resents the '2<j' error. Because of the actual shape of the 
A^-hypersurface, some of the error bars are highly asym- 
metric and the '2cr' error is often very different from twice 
the 'lcr' error. The black band at the bottom of each plot 
indicates the estimated systematic error. Overall normal- 
ization error is not included in the systematic-error band. 
These, and other results presented below are available in 
table form[24|. 

Table ITT1 shows how each of the 19 directly measurable 
spin-correlations (and Iq) can be calculated from the pa- 
rameters a for the spin-scattering matrix. While the spin- 
correlations can be calculated directly from the best-fit 
values of the parameters, as explained above, their er- 
rors cannot be found by propagation of the errors on the 
parameters. Results for I (= < > ) have already 
been shown in Fig. above. Results for the polarization 
Q[^a] = Q[ n ^\^ often denoted Pa and P A , are shown in 
Fig. El as filled circles with 'lcr' and '2cr' errors indicated. 
Also shown are results from the previous PS 185 measure- 
ment Q at pp — 1.642 GeV/c which can be seen to be 
in good agreement. Similarly Fig. HUIshows those 2-spin 
correlations which can be measured without a polarized 
target. These correlations of the spins of the final-state 
A and A are commonly denoted C mmi C m i, C nn , and C'u. 
Again good agreement is seen with the previous results. 

A particular combination of these observables, which 
has direct physical interpretation, is the singlet fraction 
(the fraction of AA pairs which are produced in a spin- 
singlet state) which can be written as 

S F = j(l - <2[n A ,n A ] + Q[ra A ,m A ] + Q[Ia,1\]) (13) 

This can be calculated directly from the spin-scattering 
matrix parameters as 



This is shown in Fig. ^2 Again the errors have been as- 
signed by directly determining the limits of change in Sf 
for a maximum acceptable change in log-likelihood. Pre- 
vious PS185 results |2| at pp = 1.642 GeV/c are shown 
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FIG. 7: Fit results for spin matrix parameters a, b, and c. The arbitrary phase is chosen by constraining the parameter a 
to be real and non- negative. Statistical errors are shown on each data point, with 2<r error bars superimposed (dashed). The 
estimated systematic error width is shown at the bottom of each plot (dark-shaded region). 



for comparison. The earlier results were determined from 
Eq. (| 1 31) using spin-correlations which had been sepa- 
rately extracted from the data. Unphysical negative val- 
ues could then occur as a result of statistical fluctuations 
or heightened sensitivity to systematic errors in the lin- 
ear combination of observables. The present results were 
extracted using Eq. (|14|) and so are constrained to be 
non-negative throughout the range of their error bars. It 
is interesting to note that the often-accepted empirical 
rule that Sf = for this reaction is clearly broken at 
back angles. 

Figure IT2l shows results for Q[n p ,n\\ and Q[n p ,n]^\ 
conventionally denoted as D nn and K nn , respectively. 
These results, which have already been published in [14j . 
are seen to disagree strongly with predictions from both 
meson-exchange (?) and quark-gluon models. While 
these results have been published, the present paper is 
the first to document the details of the technique used 
to extract them. Measurement of these spin-correlations 
was the main goal of this experiment because two com- 
peting classes of models made differing firm predictions 
for these observables. While both classes of model 
had enjoyed success in explaining the observations made 
with unpolarized targets, these results suggest that ad- 
ditional dynamics will have to be included into the mod- 
els. The wealth of additional spin-dynamics information 



presented below may help constrain and test refinements 
made to match the surprising results in these two spin- 
correlations. 

The remaining 12 directly measurable spin-correlations 
are shown in Fig. The first of these, Q[n p ], is the 
analyzing power, often denoted A„. The remainder are 
correlations between initial proton spin and components 
of one or both final-state spins. Errors are seen to be 
small enough, even on most 3-spin correlations, to allow 
structures to be clearly resolved. This underscores the 
advantage of fitting the spin-scattering matrix parame- 
ters. If these directly measurable spin-correlations had 
been determined by a direct fit of Eq. (|11|) to the ob- 
served distribution, their errors would have been much 
larger and so meaningful structure would have been im- 
possible to extract in most cases. 

Of the 256 spin-correlations, Q\jp, k p , fi^, va], defined 
by Eq. 10, one is trivially unity and 128 are constrained 
to be zero by parity conservation of the strong interac- 
tion. An additional 88 can be neglected because symme- 
try requires that they are identical to (or the negative 
of) another one which is being considered. In this sense 
there are 39 non-trivial spin correlations in addition to 
Iq. These cannot be said to be forty independent ob- 
servables since they can be expressed in terms of just 
eleven real parameters of the spin-scattering matrix at 
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FIG. 8: Fit results for spin matrix parameters d, e, and g. Statistical errors are shown on each data point, with 2a error bars 
superimposed (dashed). The estimated systematic error width is shown at the bottom of each plot (dark-shaded region). 



measurable in the present experiment and have been pre- 
sented above. However, since the spin-scattering matrix 
is fully determined, the remaining 20 spin-correlations 
can equally well be extracted just as the directly mea- 
surable ones are in the present analysis. The eight such 
observables shown in Fig. [21 would be directly measur- 
able (i.e. would appear explicitly in the description of the 
angular distribution) without a polarized beam if the tar- 
get were longitudinally polarized. Although the target 
polarization in the present experiment is purely trans- 
verse, these observables are still determined, and in some 
cases determined quite accurately, in this present mea- 
surement. Similarly, Fig. 1151 gives the results for spin ob- 
servables which would be directly measurable only if the 
antiproton beam were polarized. Again, some of these 
are quite well determined by the present data set. In 
rare cases a double-minimum in the log-likelihood func- 
tion, Ai, results in disjoint regions falling within the la 
limit. These are indicated in Fig. 1151 bv a second disjoint 
error bar. 



VII. DISCUSSION AND CONCLUSIONS 
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FIG. 9: Current results for polarization, Q[tia] = Q[fi5.] 
(filled circles). The previous measurement 0] at 1.642 GeV/c 
is superimposed (open squares). Statistical errors are shown 
on each current data point, with 2a error bars superimposed 
(dashed). The estimated systematic error width is shown at 
the bottom of each plot (dark-shaded region). 



each c . m .. However there are 40 quantities which would 
be directly measurable given arbitrary beam and target 
polarization. Of these, twenty (including Iq) are directly 



The method of determining the spin-scattering matrix, 
suggested in ^1 uas been successfully applied in practice. 
This is a unique case in which the full spin structure of a 
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FIG. 10: Current results for previously measured spin correlations between A and A (filled circles). Previous measurements 
are superimposed (open squares). Statistical errors are shown on each current data point, with 2a error bars superimposed 
(dashed). The estimated systematic error width is shown at the bottom of each plot (dark-shaded region). The spin correlations 
shown are a) Q[m A ,mA], b) Q[mx, Ja], c) Q[ti a ,iia], and d) Q[I a ,Ia\- 
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FIG. 11: Current results for singlet fraction Sf (filled cir- 
cles). The previous measurement |2[ is superimposed (open 
squares). Statistical errors are shown on each current data 
point, with 2a error bars superimposed (dashed). The esti- 
mated systematic error width is shown at the bottom of each 
plot (dark-shaded region). 



two-fermion interaction has been determined from a sin- 
gle measurement. The self-analyzing property of hyper- 
ons combined with a transversely polarized target allows 
this unusual access to the spin structure of the produc- 
tion of strange-anti-strange quark pairs. The data set of 
about 2000 events per O c m .-bin has proven sufficient to 
accurately determine the parameters and from them to 
learn the spin-correlations as well as other functions such 
as Sf and < |^ >. By construction, the results are guar- 
anteed to be internally consistent, obeying all constraints 
imposed by the symmetries of the strong interaction un- 
der parity and charge conjugation. Numerical values for 



these results are available |24j . 

Apart from the small background subtraction at back- 
angles, all results shown for each angular bin have been 
obtained completely independently of the results at other 
angular bins and are based on non-overlapping sets of 
events. The smooth variations as a function of c . m . seen 
in most of the spin-correlations is in no way built-in to the 
analysis technique. The fact that the angular variation 
of the data appears smooth is a reassurance that the 
entire chain of event reconstruction and data analysis is 
performing reasonably and extracting meaningful results. 
Similarly the fact that the bin-to-bin 'scatter' in the data 
appears to be consistent with the assigned error bars is 
an independent verification of the validity of the method 
of error analysis. 

The most important aspects of these results, relating 
to Q\n„, n\] and Q[n p ,n A ] have already been presented 
in [lj]. As shown in Fig. El the measured values of 
these spin-correlations differ markedly from the predic- 
tions of a meson exchange (MEX) modelp} (solid line) 
and a quark-gluon-inspired (QG) model [£j (dotted line) 
despite the fact that both these models reasonably de- 
scribe the significant spin-structure observable with an 
unpolarized target. 

All MEX models generally predict a large tensor inter- 
action which couples spin-triplet pp initial state to spin- 
triplet AA final state, flipping the spin in the process. For 
this reason, both the spin transfer, Q[n p , n A ], and the de- 
polarization, Q[n p , tia], are predicted to be strongly nega- 
tive. This prediction holds even in the presence of initial- 
and final-state interactions, which have been included in 
the prediction shown by solid lines in the figures. The 
measured values are far less strongly negative than the 
predictions and, in fact, are positive at forward angle in- 
dicating that both final-state particles tend to be aligned 
with the initial proton spin. Furthermore, Q[n p ,n A ] is 
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FIG. 12: Results for spin transfer observables Q[n p , tia] (often called depolarization, D nn ) and Q[n p ,riA] (often called spin- 
transfer, Knn) at 1.637 GeV/c, compared to MEX model prediction (solid) and QG model prediction |?| (dashed) at 
1.642 GeV/c. Statistical error bars are shown on each data point, with 2a error bars superimposed (dashed). The estimated 
systematic error width of the measurement is shown at the bottom of each plot (dark shaded region). 



strongly positive at back angles meaning that the nor- 
mal component of the proton's spin is transferred to the 
A in contrast to the MEX prediction. 

All existing calculations using QG models have been 
restricted to 3 Po, 'vacuum' terms, and 3 <Si 'gluon' terms. 
So the interaction is purely spin triplet, having Sf = 
built into the model. Here there is a much smaller tensor 
interaction and so less spin-flip. It was because of this 
characteristic difference between QG and MEX models 
that Q[n p ,nA] was first suggested [ll| as an observable 
which would distinguish experimentally between the two 
classes of model. A vanishing singlet fraction necessarily 
implies that Q[n p ,nA] = Q[n p ,riA}- This is reflected in 
the predictions of the QG model shown as dotted lines in 
Fig-El As shown in Fig. ^2 however, the singlet fraction 
is distinctly non-zero at back angles, in violation of the 
assumptions of the existing models. This manifests itself 
in the data as a large difference in back-angle behavior 
of the two distributions. 

Apart from large 6 c . m . and very small c .m., the QG 
model is seen to dramatically over-predict both Q [n^, ua] 
and Q[n p , n A ]- Predictions have also been made |25j| for 
the transfer of other components of proton spin to the 
spin of the A and A. These are compared to the data 
in Fig. El While the disagreement is not as striking, 
partly because of the relatively larger error bars, it is 
clear that significant modification of the model will be 
needed to match these correlations and to predict the 



others reported here. A modification which is clearly 
required is inclusion of singlet strength, such as 1 So. 

While a great wealth of information has been gained 
on the spin dynamics of pp — > AA at p p = 1.637 GeV/c, 
few conclusions can be drawn because theoretical models 
lag significantly behind in understanding the data at this 
point. Availability of this data may inspire increased 
theoretical activity. 

Among the most precisely determined spin correlations 
is one which has not previously been measured. The an- 
alyzing power (a correlation of spin only with scattering 
angle, not with other spins) Q[n p ] (Fig. IT5fal. usually 
denoted A n , is constrained to vanish at cos(0 c m .) = 
±1 but has now been determined to be strongly posi- 
tive in the forward hemisphere and mostly negative for 
back angles. The complex angular structures of the 2- 
spin correlations, Q[toa, toa], Q[m A ,lA], Q["-A! n A], and 
Q[/a,^a], (Fig. El a through d, respectively) although 
consistent with previous measurements Q, is now re- 
vealed in far more detail. Among the 3-spin correlations, 
some of the directly measurable ones show the cleanest 
structure with Q[n p , toa, ^a] (Fig. 113k) and Q[n p , I a, toa] 
(Fig. I13L ) showing remarkably similar behavior, while 
it is not trivially expected that they should be equal. 
A distinctly different structure is seen in Q[n p , toa, toa] 
(= — Q[n pi Ia 1 I a]) (Fig- El) which is more compressed 
to forward angles. Even some of the 4-spin correla- 
tions are determined well enough to unveil distinct an- 
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FIG. 13: Results for the twelve additional spin observables which appear directly in the measured angular distribution. 
Statistical and systematic error estimates are displayed, as above. The spin observables displayed are a) Q[n p ], b) Q[m p ,mA], 
c) Q[m p ,lA], d) Q[m p ,m A ], e) Q[m p ,l A ], f) Q[m p ,m A , tia], g) Q[m p ,n A ,m A ], h) Q[m p ,n A ,l A ], i) Q[m p , l A , ua], j) Q[n p , m A , 771a], 
k) Q[n p ,m A ,l A ], and 1) Q[np,/ A ,m A ] 



gular structure. Regions of non- vanishing correlation are 
seen for example in Q[mp, l p , m A , m\], Q[m,p, l p , m A , l A ], 
and Q[mp, l p , l A , m\\ (Figs. IT5h. II 51. and 115k. respec- 
tively). It would be fruitless to speculate on the meaning 
of each of these structures individually. A coherent pic- 
ture will required theoretical modeling to simultaneously 
explain all available spin-correlations or equivalently to 
directly predict the coefficients of the spin-scattering ma- 
trix. This difficult task is made all the more difficult 
by the strength of the initial- and especially final-state 
interactions. But even in their presence, the physics re- 
duces to just eleven real parameters at each angle. There 



may be advantages to comparing theories to the experi- 
mentally determined spin-matrix parameters. They may 
tie in more directly to the underlying spin physics of 
the model. Also the process of determining the quality 
of the agreement is simplified by removing the 'double- 
counting' which is inherent in comparing up to 40 spectra 
when all the physics reduces to just eleven sets of param- 
eters. 
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FIG. 14: Results for eight spin observables which do not appear directly in the measured angular distribution. Direct 
measurement of an individual one of these observables would require a longitudinally polarized target proton but would not 
require beam polarization. Statistical and systematic error estimates are displayed, as above. The spin observables displayed 
are a) Q[l p ,m A ], b) Q[l p ,l A ], c) Q[l p ,m A ], d) Q[l p , m A , n A ], e) Q[l p , n A , m\], f) Q[l p ,n A ,l A ], g) Q[l P ,l A ], and h) Q[l p ,l A ,n A ]. 
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